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Abstract 

Background: A number of cellular- and molecular-level studies of autophagy assessment 
have been carried out with the help of various biochemical and morphological indices. 
Still there exists ambiguity for the assessment of the autophagy status and of the causal 
relationship between autophagy and related cellular changes. To circumvent such 
difficulties, we probe new quantitative indices of autophagy which are important for 
defining autophagy activation and further assessing its roles associated with 
different physiopathological states. 

Methods: Our approach is based on the minimal autophagy model that allows us 
to understand underlying dynamics of autophagy from biological experiments. 
Specifically, based on the model, we reconstruct the experimental context-specific 
autophagy profiles from the target autophagy system, and two quantitative indices 
are defined from the model-driven profiles. The indices are then applied to the 
simulation-based analysis, for the specific and quantitative interpretation of the system. 

Results: Two quantitative indices measuring autophagy activities in the induction of 
sequestration fluxes and in the selective degradation are proposed, based on the 
model-driven autophagy profiles such as the time evolution of autophagy fluxes, levels 
of autophagosomes/autolysosomes, and corresponding cellular changes. Further, with 
the help of the indices, those biological experiments of the target autophagy system 
have been successfully analyzed, implying that the indices are useful not only for 
defining autophagy activation but also for assessing its role in a specific and 
quantitative manner. 

Conclusions: Such quantitative autophagy indices in conjunction with the computer- 
aided analysis should provide new opportunities to characterize the causal relationship 
between autophagy activity and the corresponding cellular change, based on the 
system-level understanding of the autophagic process at good time resolution, 
complementing the current in vivo and in vitro assays. 

Keywords: Autophagy, Quantitative indices, Minimal autophagy model, 
Computer simulations 



Background 

Macroautophagy (hereafter referred to as autophagy) is a key homeostatic mechanism 
for the turnover of such intracellular components as proteins/organelles [1], and is fur- 
ther related to various human diseases such as cancer, metabolic disorders, and neuro- 
degenerative diseases [2-7]. This has brought forth a number of physiological and 
molecular-level studies of autophagy in the last decades. 
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In the process of autophagy, abnormal and/or resident proteins/organelles degrade 
into representative metabolic/energy precursor molecules, amino acids and ATP, which 
can be used as new building blocks and energy sources, respectively. Specifically, the 
autophagic process begins with the formation of initial double-membrane structures 
called autophagosomes containing sequestered proteins/organelles. Then the autopha- 
gosomes fuse with endosomes/lysosomes to form autolysosomes. Finally, the contents 
engulfed in the autolysosomes are hydrolyzed via intralysosomal hydrolysis. 

The operation mechanism of the process is extremely complex. The three consecu- 
tive steps, i.e., autophagosome formation, autolysosome formation, and intralysosomal 
hydrolysis, are operated independently, exhibiting qualitatively/quantitatively different 
responses to different intra- or extra-cellular perturbations. Hence, assessment of the 
autophagy activity via specific markers such as autophagosomes and autolysosomes 
often leads to under-/over-estimation of the autophagy activity. For example, as re- 
vealed in experiment [8], if the rate of autolysosome formation exceeds that of autopha- 
gosome formation, the steady-state concentration of autophagosomes could be 
detected less, leading to an underestimation of the autophagy status. Moreover, the au- 
tophagic process is regulated by the complex positive- negative feedback mechanism 
(see Figure 1): The output of the process, i.e., recycled amino acids and ATP, acts as 
the input of the system, which independently manipulates each of the steps in a 
concentration-dependent manner [9-14]. Such feedback loops may be beneficial for the 
delicate balance of the cellular homeostasis. In the perspective of causality analysis/ 
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Figure 1 Schematic representation of the model system and quantitative indices. The model 
assumes a three-compartment description of the autophagic process: protein/organelle, autophagosome, 
and autolysosome compartments. Solid and dotted arrows denote the autophagic feedback loops among 
the compartments and non-autophagic degradation, respectively. The autophagy-related rates, including 
protein/organelle synthesis R s {C Al C 0 ), autophagosome formation R g] {C Al C a ) and R g2 {C Al Ca), autolysosome 
formation intralysosomal hydrolysis Rh(C A ), and non-autophagic degradation R d] and R d2 are described 
in Appendix A. Dashed arrows indicate the quantitative indices f ci O s , and 0, of protein/organelle quality, 
autophagy selectivity, and autophagy induction, respectively. 
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assessment, however, it could lead to difficulties in interpreting the cause-and-effect re- 
lationship between the autophagy activity and the concomitant cellular change. 

Even if those uncertainties in the assessment of the status and causality of autophagy 
are removed, there still remains ambiguity as to the roles of autophagy in the human 
disease since autophagy exhibits dual effects on the development and progression of 
various human diseases in a context- and activation-degree-dependent manner. Espe- 
cially, autophagy plays a dual role in the tumor cell viability [15]: In some cases autoph- 
agy prevents or suppresses tumor progression whereas in other cases autophagy can 
also accelerate tumorigenesis or promote the survival of tumor cells. Under energy- 
deficient conditions, autophagy is usually activated for the rapid supply of essential en- 
ergy/metabolites to promote cell survival. In contrast, it is not rare that induced au- 
tophagy contributes to apoptosis/necrosis under various cellular perturbations other 
than the energy deficiency. Accordingly, it is intriguing to interpret the newly eluci- 
dated mechanism of autophagy: whether the induction of autophagy contributes to the 
prevention of the disease and how it works. In particular, evolutionally new molecular- 
level studies sometimes contradict existing studies. 

Based on the recent molecular- and physiological-level studies of the mechanism and 
the role of autophagy in the human disease, on the other hand, there is research being 
carried out on the development of the treatment methods or drugs that can potentially 
regulate or control autophagy [16-19]. In view of these, it is very desirable to develop 
specific and quantitative indices based on adequate mathematical model which could 
provide a set of reliable criteria for defining the autophagy status and further the kinet- 
ics of the process. To date, various indices based on biochemical [20-25] and morpho- 
logical methods [26-30] for the detection of the autophagic sequestration of proteins/ 
organelles and for the measurement of the turnover of autophagic compartments and/ 
or autophagy-related markers have been developed for the assessment of autophagy. 
However, those indices are not quite satisfactory for the specific and quantitative as- 
sessment of the autophagy pathway as well as the system-level interpretation of the 
process, especially in mammalian cells. 

The aim of the present study is to propose the experimental context-specific quan- 
titative indices of the autophagy activity based on the minimal autophagy model and 
to apply the indices to the simulation-based analysis of the target autophagy system 
(see Figure 2). Specifically, based on the minimal autophagy model [31], we recon- 
struct underlying profiles of the autophagic process such as the time evolution of au- 
tophagy fluxes, levels of autophagosomes/autolysosomes, and corresponding cellular 
changes at good time resolution, from the biological experiments of the target au- 
tophagy system [12-14,26,32,33]. Then, two quantitative indices measuring autophagy 
activities in the induction of sequestration fluxes and in the selective degradation are 
proposed from the model-driven profiles. Finally, the biological experiments are 
interpreted with the help of the indices, not only to examine how the autophagy sys- 
tem responds to cellular damaging but also to prove the causal relationships among 
the steady-state autophagy levels, autophagic fluxes, and corresponding cellular 
changes in a specific and quantitative manner. Developing such quantitative autoph- 
agy indices from the minimal model will be highly challenging but indispensable for 
the assessment of autophagy activity and of its roles associated with different physio- 
pathological states. 
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Figure 2 Development of the experimental context-specific quantitative autophagy indices from 
the minimal model. The workflow of the process is as follows: 1. Establishment of the minimal autophagy 
model for the target autophagy system. 2. Reconstruction of underlying dynamics of the autophagic process at 
good time resolution. 3. Introduction of two quantitative autophagy indices measuring autophagy activity 
based on the model-driven autophagy profile. 4. Computer-aided analysis of the biological data with the help 
of the indices. 5. Specific/quantitative interpretation and system-level understanding of the system. 



Methods 

Model systems and simulations 

We first outline briefly the minimal model used in developing quantitative indices of 
the autophagy activity [31]. In this model, the multi-step autophagic process is di- 
vided into autophagosome formation, autolysosome formation, and intralysosomal 
hydrolysis steps, and the intracellular ATP and amino acids are considered to be the 
key molecules serving as a bridge among the autophagic process, non-autophagic 
degradation, and protein/organelle synthesis (see Figure 1). The dynamics of the 
model system are described by eight coupled differential equations 1 to 6, which are 
solved via the 5th order Runge-Kutta method for very high precision. Since the 
present study is designed to analyze a specific autophagy mechanism rather than to 
examine the general principle of autophagy, we have focused on the target autoph- 
agy system in the parameter selection: The key parameters used in the simulations 
have been fixed in accord with carefully selected biological data of the target au- 
tophagy system [12-14,26,32-34]. A few adjustable parameters, which have been set 
to reflect the experimental context, have minor influence on the simulations (see 
Table 1). 



Dynamic equations 

In this model, variations of the autophagosome concentration with time are determined 
by the difference between the autophagosome formation specific rate R gi and the auto- 
lysosome formation specific rate R t . Denoting by C gi the concentration of autophago- 
somes originating from proteins/organelles % we write the equations for the dynamics 
of autophagosomes in the form (i = 1, 2): 
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Table 1 Parameters in the model 



Symbol 


Definition 


Unit 


Value 


Reference 


9 


Rate constant for autophagosome formation (normal value) 


s" 1 


* 1.12 x 10" 5 


[12], [13], [14] 




Constant for autophagosome formation (ATP dependency) 


mM 


* 4.01 


[12], [14] 


Pn 


Constant for autophagosome formation (ATP dependency) 


mM 


* 3.00 


[12], [14] 


9 


Constant for autophagosome formation 
(amino acids dependency) 


mM 


4.50 




Yg 


Constant for autophagosome formation 
(amino acids dependency) 


(unitless) 


* 1.22 


[13] 




Constant for autophagosome formation 
(amino acids dependency) 


H A- 1 

mM 


* 7.49 x 1 0" 2 


[13] 


n 


Rate constant for autolysosome formation 


s- 1 


*2.47x10" 5 


[12], [14] 


ki 


Constant for autolysosome formation (ATP dependency) 


mM 


* 4.01 


[12], [14] 


Pi 


Constant for autolysosome formation (ATP dependency) 


mM 


* 3.00 


[12], [14] 


r h 


Rate constant for intralysosomal hydrolysis 


s" 1 


* 1.39X10" 5 


[12], [14] 


6 h 


Exponent for intralysosomal hydrolysis (ATP dependency) 


(unitless) 


* 7.24X10" 1 


[12], [14] 


k h 


Constant for intralysosomal hydrolysis (ATP dependency) 


mM 


* 2.99 


[12], [14] 


h 


Rate constant for protein/organelle synthesis 


mM • s" 1 


* 1.48X10" 5 


[34] 


k 5 


Constant for protein/organelle synthesis 
(amino acids dependency) 


mM 


* 1.77 x 10 1 


[34] 


Am) 


ATP concentration corresponding to maximal 
protein/organelle synthesis rate 


mM 


3.00 





Parameters with asterisks (*) are fixed, determined from biological experiments [12-14,34]. Those without asterisks are 
adjustable, depending on experimental setups. 



dCgi 
dt 



(i) 



where C S i represents the concentration of S/. Next, the intracellular concentration C# 
of autolysosomes, originating from S; is determined by the difference between the auto- 
lysosome formation specific rate Ri and the intralysosomal hydrolysis specific rate R^. 
The equations governing the dynamics thus read (i = 1, 2): 



dCu 
dt 



Ri(t-r) C gi (t-T)-R h C u 



(2) 



Note that the autolysosome concentration at time t is affected by the autophagosome 
concentration at time t - r, earlier by the delay time r which is taken to be 8 minutes 
(r = 480 s) [26,32,33]. 

We have defined resident proteins/organelles Si as the proteins and organelles which 
conduct normal functions in the cell, and assumed that they are translated from normal 
folding intermediates transcribed from DNA normally into RNA. On the other hand, 
by abnormal proteins/organelles S 2 , we have meant the proteins and organelles which 
conduct abnormal functions in the cell and assumed that they are made from two dis- 
tinct sources: either from misfolded proteins and peptides, caused by genetic variants 
and mutations or intracellular conditions, or from resident proteins and organelles, 
damaged or aged by harmful conditions. Incorporating these, we have described the dy- 
namics of Si and of S 2 by the evolution equations for the concentrations C S i and C 52 
of Si and S 2 , respectively: 
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dC S i 



(l-a)R s -R dl -pCsi-R g iC S i 



(3) 



dt 



dC S 2 



aR s - Rd2 + PC S \ - R g2 C S 2, 



(4) 



dt 



where R s represents the (total) protein/organelle synthesis rate (from DNA) and a is 
the fraction of S 2 in the protein/organelle synthesis. Accordingly, Si and S 2 are pro- 
duced at the rates of (1 - a)R s and aR s , respectively. Further, R di represents the non- 
autophagic degradation rates of S/ (for i=l,2) and ft is the specific rate of deterioration 
of Si, i.e., transformation from Si to S 2 . 

The dynamics of intracellular amino acids, the concentration of which is denoted by 
C at takes the form 



of which the first and second term on the right-hand side correspond to the increase of 
amino acids due to intralysosomal hydrolysis and nonautophagic protein/organelle deg- 
radation, respectively, with appropriate constants p a and pi d describing the mean num- 
bers of amino acids produced from autophagic degradation and from non-autophagic 
degradation, respectively. The third term represents the net intracellular amino acid 
generation rate due to various intracellular metabolisms other than autophagy, and has 
been defined according to R a = p c R s with an appropriate constant ft c . The last term stands 
for the reduction of amino acids due to protein/organelle synthesis, with the constant p s 
denoting the mean number of amino acids within a protein/organelle. 

It has been assumed that the intracellular ATP concentration C A increases due to 
intralysosomal hydrolysis and non-autophagic protein/organelle degradation via cytosolic 
and mitochondrial ATP production. The corresponding dynamic equation for ATP reads: 



where v a and v d describe the mean numbers of ATP molecules produced from autoph- 
agic degradation and from non-autophagic degradation, respectively, and v s gives the 
mean number of ATP molecules consumed in unit protein/organelle synthesis. The net 
intracellular ATP generation rate R A is given by the difference between the cytosolic 
and mitochondrial ATP production rate and the ATP consumption rate, and is as- 
sumed to be R A - v c R s with a constant v c . 

For simplicity, we have supposed that an average protein/organelle in the virtual cel- 
lular system is composed of 500 amino acid residues, implying that 500 amino acids 
are consumed for the synthesis of a protein/organelle. Further, considering the fact that 
elongation of one amino acid during translation requires approximately four ATP mol- 
ecules, we have assumed 2000 molecules of ATP involved in the synthesis of a protein/ 
organelle. However, since the efficacy of protein/organelle recycling is expected to be 
less than 100%, the increases of amino acids and ATP due to intralysosomal hydrolysis 
or non-autophagic protein/organelle degradation should be less than 500 and 2000 
molecules per degradation of one protein/organelle, respectively. In this study, we have 




(5) 




(6) 
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chosen the base run parameter values as pi a = \i d = v a = v d = 300, \i c = 200, v c = 1700, \i s = 
500, and v s = 2000. 



Results 

Quantitative indices from minimal models 

Solving Eqs. 3 and 4, we have obtained the time evolution of concentrations C 51 and 
C S 2 of resident Si and abnormal proteins/organelles S 2 , respectively (see Figure 3). They 
display oscillations with the period of 12 h 55 min (or the natural frequency 0.0774 h" 1 ) 
[35-39], with Csi and Cs2 of the average values of C$i and C 52 given by 8.57 mM and 
2.52 mM, respectively. 

Based on such model-driven profiles of intracellular concentrations of proteins/ or- 
ganelles, we propose the fractional abnormal protein/organelle concentration / c as a 
simple index of the cellular protein/organelle quality [40]: 



Cs2 



Csi + C 



52 



(7) 



which takes values between 0 and 1: While the value f c ~ 0 addresses that most of the 
proteins/organelles in the system are normal, f c « 1 indicates that the majority of the 
proteins/organelles are damaged. 
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Figure 3 Time evolution of the protein/organelle concentrations. Red solid and blue dotted lines 
describe the time evolution of the concentrations C 51 (top) and C 52 (bottom) of resident and abnormal 
proteins/organelles, respectively. Data have been obtained at r g = rf' and /3 = /3 (0) . 
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Figure 4 exhibits the time evolution of concentrations C g \ and Cgz of autophago- 
somes from Si and from S 2 > respectively (left axis), and the time evolution of concen- 
trations C/i and C/2 of autolysosomes from Si and from S 2 , respectively (right axis). 

It has turned out that they display oscillations with the period of 12 h 55 min (or the 
natural frequency 0.0774 h" 1 ) in the absence of external driving [41-48]. The percentile 
ratios of the peak-to-peak amplitudes to the mean levels of the oscillations are 0.929% 
for C g i, 0.0676% for Cg2, and 20.8% for Qi and for Q 2 , with negligibly small standard 
deviations (of the order 10~ 5 relative to the amplitudes). Note here that the oscillation 
amplitudes of autophagosome concentrations are far smaller than those of autolyso- 
some concentrations. The average concentrations C gX and Cgi turned out to be 
0.0589 mM and 1.90 mM, respectively and those of autolysosomes, Qi and Q 2 , to be 
0.0156 mM and 0.503 mM, respectively. 

Furthermore, the minimal autophagy model allows us to characterize quantitatively 
the relationship between the autophagy activity and the cellular change, based on the 
system-level understanding of the autophagic process at good time resolution, which 
may not be possible via the current in vivo and in vitro measurements. For example, 
from the model-driven profiles of the steady-state concentrations of autophagic inter- 
mediates (Figure 4), together with the dynamics of proteins/organelles (Figure 3), 
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Figure 4 Time evolution of the autophagosome and autolysosome concentrations. Red solid and 
blue dotted lines describe the time evolution of the autophagosome concentrations C g] and C g2 and 
autolysosome concentrations C n and C /2 , respectively, transformed from S] (top panel) and from S 2 



(bottom panel). Data have been obtained under the condition r g = ry and /3 = 0.2 %/h = 0 
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substrate selectivity of the autophagic process has been observed. There are significant 
differences between the production rates of autophagosomes/autolysosomes from resi- 
dent proteins/organelles and those from abnormal proteins/organelles: The steady-state 
concentration ratio of Qi to C 51 is given by 0.00182 whereas that of Q 2 to C S2 is 0.200. 
In addition, the ratios of C gX to C 51 and of Cg 2 to C S2 read 0.00687 and 0.754, respectively. 
For a more specific and quantitative analysis of various facets of the substrate selectivity 
of autophagy [40,49-54], we here propose a quantitative index of the autophagy selectivity 
(see Figure 1). 

The autophagy selectivity index O s , which is designed to quantify the selective au- 
tophagic degradation of abnormal or resident protein/organelle, is given by: 

Os^-Bl, (8) 

where Cu denotes the average concentration of autolysosomes from S/(z = 1, 2). Note 
that positive values of O s imply that autophagic degradation of S 2 exceeds that of Si 
whereas negative values correspond to the degradation of Si exceeding that of S 2 . 

This model allows us to reconstruct not only the dynamics of the steady-state con- 
centrations of autophagosomes/autolysosomes but also the time evolution of autophagy 
fluxes. As shown in Figure 5, fluxes of sequestration, of maturation, and of intralysoso- 
mal hydrolysis from S 2 , denoted by F s2 , F m2 , and F h2 , respectively, and those from Si, 
denoted by F sl , F ml , and F hv display synchronized oscillations. Note also that the values 
of fluxes from S 2 are much greater than those from Si. 

Based on such comprehensive information as to the kinetics of the process obtained 
from the minimal model, we propose the autophagy induction index O x , which is de- 
vised to measure the induction of the total sequestration flux coming from both resi- 
dent proteins/organelles Si and abnormal proteins/organelles S 2 : 

Oi= F sl +F s2 , (9) 
where F$i denotes the average sequestration flux for = 1, 2). 



Quantitative analysis through the use of autophagy indices 

To determine whether the indices given in the previous section are useful for the spe- 
cific and quantitative assessment of the autophagy pathway and for the system-level un- 
derstanding of the process, we in this section apply the indices to the interpretation of 
the target autophagy system, metabolically controlled autophagic protein/organelle deg- 
radation in the rat hepatocyte [12-14,26,32,33]. To be specific, we conduct the 
simulation-based analysis of how the autophagic process responds to cellular damaging 
and prove the causal relationships among the steady-state autophagy levels, autophagic 
fluxes, and corresponding cellular changes with the help of the indices. 

It is shown in Figure 6 that O x grows drastically with the (specific) deterioration rate 
ft until the value ft ~ 0.12, which illustrates autophagy induced against the cellular dam- 
aging rate. As is raised beyond the value /?- 0.12, O x increases gradually and displays 
a plateau. On the other hand, as is increased from zero, O s stays at relatively high 
positive values until ~ 0.12, beyond which it reduces gradually. 

Therefore, both the non-selective mode of sequestration fluxes (represented by <£>{) 
and the selective mode of the autophagic degradation of abnormal proteins/organelles 
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t (hour) 

Figure 5 Autophagic fluxes in sequestration, maturation, and intralysosomal hydrolysis. Red solid, 
blue dotted, green dashed lines describe the time evolution of sequestration fluxes F s1 and F s2l maturation 
fluxes F m] and F m2l and intralysosomal hydrolysis fluxes F h] and F h2l respectively, originating from S] (top) 
and from S 2 (bottom). Data have been obtained at r g = and /3 = /3 (0) . 



(described by O s ) have been evaluated in a specific and quantitative manner with the 
help of the indices: The increasing behavior of <&i with the deterioration rate ft suggests 
that the total sequestration flux coming from both resident and abnormal proteins/or- 
ganelles is induced, resisting against the cellular damaging level In addition, the posi- 
tive values of O s in the range of ft from 0 to 0.5 (%/h) indicates that abnormal 
proteins/organelles are selectively removed via autophagy. 

We then analyze how the promotion or suppression of autophagy affects the cellular 
quality control. Figure 7 exhibits the fractional abnormal protein/organelle concentra- 
tion f c depending on the autophagy indices O t and O s , in response to varying the rate 
constant r g for autophagosome formation in Eqs. Al and A2 in Appendix A. Data have 
been obtained at the specific deterioration rate ft = 0.2(%/h) = /? (0) , where the resident 
protein/organelle synthesis rate is approximately equal to the abnormal one. 

In the case of no autophagic flux in the system (Oj = 0), the fractional concentration 
is positioned at a very high level {(f c ~ 1), indicating that most of the proteins/organelles 
in the system are damaged (top panel). As autophagosome formation is promoted, 
namely, as r g is increased, however, <&i begins to increase. In particular, the abrupt increase 

in ®! at r g = 0.03, together with the drop in f c at the same rate constant, appears 
here as the apparently discontinuous change in f c at the values between <&i ~ 0 and O t « 2. 
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Figure 6 Autophagy induction index O, (red pluses) and selectivity index O s (blue asterisks) versus 
the specific deterioration rate j8. Data have been obtained at r g = rf^ with /3 varied up to 0.5 (%/h). 



As ®! is increased further, f c reduces to the normal level and drops eventually close to 
zero. 

In addition, f c varies also inversely proportional to O s . In case that the autophagic 
degradation of resident protein/organelle is larger than that of abnormal protein/organelle 
(® s < 0),f c stays at an abnormally high level (f c « 1), namely, most of the proteins/organ- 
elles in the system are damaged (bottom panel). As O s is increased slightly above zero, 
there arises a discontinuous drop via which f c becomes restored quickly to the normal 
level. As O s is raised further and approaches unity,/, keeps decreasing toward zero. 

Discussion 

The minimal model for autophagy was originally developed to examine the dynamics 
of the autophagic process, describing specifically the rates at which autophagosome 
and autolysosome concentrations vary with time [31]. Based on the model, we have re- 
constructed underlying profiles of the autophagy process such as the time evolution of 
autophagy fluxes, levels of autophagosomes/autolysosomes, and corresponding cellular 
changes from the target autophagy system, in which the corresponding biological ex- 
periments [12-14,26,32,33] provide information only as to the changes before and after 
certain experimental perturbations on specific processes. We have then proposed quan- 
titative indices of the autophagic process, and used them to analyze how the autophagic 
degradation compensates cellular damaging. Emphasis has been paid on the stressful 
conditions, specifically, at extremely high rates of protein/organelle deterioration. As 
discussed in Figures 6 and 7, indices O x and O s have successfully characterized the 
non-selective induction rate of autophagy and selective intralysosomal hydrolysis, respect- 
ively, in the presence of physiological perturbations such as variations of the cellular dam- 
aging rate and promotion or suppression of autophagosome formation. Furthermore, as 
shown in Figure 7, f c reduces sensitively as <&i or O s is increased; this might be inferred 
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Figure 7 Fractional protein/organelle concentration f c versus autophagy indices 0| and O s . Data 
plotted with filled purple rhombi (top panel) and open blue rhombi (bottom panel) have been obtained 
via suppressing or promoting autophagosome formation at i 6 = /3 (0) . The rate constant has been raised from 
r Q = 0 to r Q = 1 0 at the increment of 0.01 (in units of ri 0 ^). 




from the result that both the induction of the total sequestration flux, represented by O x , 
and the improvement of substrate selectivity in the autophagic degradation, described by 
® s , are beneficial for the control and regulation of the cellular protein/organelle quality, 
measured by f c . On the other hand, these results may disclose the roles of evolutionary- 
conserved basal autophagy in cell survival Under extremely low levels of autophagic flux, 
the system may lose its ability to control the cellular protein/organelle quality, eventually 

resulting in cell death: As r g is reduced below r^~0.03, the fractional concentration f c re- 
mains at an abnormally high level, which may cause dysfunctions of the cell (Figure 7). 
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Particularly, it is expected that once r g reduces below rf\ normal cellular functions might 
not be recovered; this appears to be supported by the accompanying behaviors of resident 
and abnormal proteins/organelles as well as of autolysosomes, amino acids, and ATP 
(data not shown). Furthermore, in the case of no autophagic flux in the system as illus- 
trated for r g = 0, most of the proteins/organelles in the system are damaged (f c ~ 1), with 
the corresponding indices given by <&i = 0 and O s < 0 (see Figure 7). 

As the practical applications of this theoretical study, we remark implications of the 
quantitative autophagy indices for the development of a novel strategy for the assess- 
ment of autophagy. Several biochemical assays such as the measurement of autophagic 
sequestration [22], detection of the turnover of specific autophagic compartments or 
autophagy- related markers [23-25], and quantification of the autophagic protein/organ- 
elle degradation [55] have recently been suggested to provide indirect correlative data 
relating to the autophagic protein/organelle quality control. However, there still lacks 
full biochemical or molecular understanding of what distinguishes the selective and 
nonselective autophagic modes and how context- and activation-degree-dependency of 
selective/nonselective autophagy contributes to the protein/organelle quality control. 
Further, there are no absolute criteria, applicable to diverse situations, for determining 
the autophagic flux, mainly because some techniques and chemicals in certain assays 
are inappropriate, problematic, or may not work at all in other assays. Namely, it is not 
yet known whether the measurement of LC3-IIs/LC3-I and p62/SQSTMl is generally 
applicable to other cell types, for changes in LC3-II or p62 amounts are tissue- and cell 
context-dependent; this constitutes the main caveat regarding the method. Also, the 
GFP-LC3 processing assay appears to depend on cell types and culture conditions, 
which is the main limitation. Even in the same assay, when using those techniques and 
chemicals, it is critical to consider the tissue- and cell context-dependent effects. In- 
deed, in some cases, the amounts of those indirect markers may not correlate well with 
the autophagosome/autolysosome accumulation detected by electron microscopy which 
is the most reliable criterion for autophagy activity. 

In these circumstances, it is very desirable to have quantitative indices together with 
the appropriate mathematical model, which make it possible to provide a set of reliable 
criteria for the definition of the autophagy activation and further the assessment of its 
roles associated with different physiopathological states. Furthermore, the quantitative 
indices of the autophagy activation could give kinetic information as to the autophagic 
process, i.e., autophagosome formation, autolysosome formation, and intralysosomal 
hydrolysis. With such parameters available quantitatively, the worth of those conven- 
tional approaches to assessing autophagic fluxes or steady-state quantities of autopha- 
gosomes or autolysosomes could be greatly enhanced [27,56-61]. 

It is also to be noted that the proposed indices O x and O s are designed based on the 
selective profiles of autophagy, i.e., autophagic intermediates/fluxes from Si and those 
from S 2 . Further, the indices can be modified depending on the experimental settings, 
with which selective autophagy is associated [40,49-54]: In fact, Si and S 2 in this study 
can denote different proteins/organelles in specific diseases such as aging', 'intracellular 
quality control and housekeeping', and 'host defense against intracellular pathogens'. 
Accordingly, the indices should be easily applicable to various biochemical [20-25] and 
morphological [26-30] experiments, which selectively detect and quantify autophagy ac- 
tivity depending on their target substrates. Although such selective measurement of the 
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autophagy activity with high specificity has not yet been carried out within current bio- 
assay techniques [27,56,57], the indices should be useful in getting a better grasp of the 
substrate-selective autophagy activity and its role, complementing current biological 
techniques [20-30,55-57], 

Conclusion 

We have proposed quantitative autophagy indices, based on the mathematical model, 
to define autophagy activity and further assess its role. With the help of the indices in 
conjunction with the computer-aided analysis, we have characterized quantitatively the 
cause-and-effect relationship among the steady-state autophagy levels, autophagic 
fluxes, and corresponding cellular changes in response to various physiological pertur- 
bations, which may not be probed via biological assays. 

Our study acts as a natural link between experimental and computational/theoretical 
biology. Such an integrative approach should possibly lead to a comprehensive under- 
standing of the control and/or regulatory mechanism of autophagy and reduce ambigu- 
ity as to causes and effects in the complex autophagy system. From a medical 
viewpoint, it should provide us new insight into the role of autophagy in various human 
diseases, including cancer, metabolic disorders, and neurodegenerative diseases and fur- 
ther help to develop new drugs or treatment methods which target specific autophagy 
pathways [16-19]. We hope to develop more realistic model-based quantitative indices 
of autophagy activity and new methods of monitoring autophagy, complementing re- 
cent biochemical assays; this is left for further study. 

Appendix A 

Autophagosome formation 

Under normal conditions, it appears that abnormal proteins/organelles degrade preferen- 
tially via basal autophagy [49-52]. As the intracellular energy/nutrient reduces due to, e.g., 
starvation or increased metabolic demand, however, both resident and abnormal pro- 
teins/organelles are degraded non-selectively by bulk autophagy for the rapid supply of 
the essential energy/metabolite molecules. Therefore, it is assumed in this model that the 
autophagosome formation rate from resident proteins/organelles, which is lower than that 
from abnormal proteins/organelles under normal conditions, becomes gradually equal to 
that from abnormal proteins/organelles as the amino acid concentration is decreased 
[5,11,62,63]. Note, however, that the molecular mechanism of these steps is not included 
since the mechanism, via which ATP and amino acids control autophagosome formation, 
is relatively poorly understood. For example, preliminary studies show merely that amino 
acids regulate the LC3 level at the mRNA expression level [64] . Instead, we have obtained 
the dependence of these steps on intracellular ATP and amino acid concentrations, based 
on physiological-level experiments [12-14] (see Figure 1). 

In consequence, we thus take the autophagosome formation specific rates R g i (from resi- 
dent proteins/organelles Si) and R^ (from abnormal proteins/organelles S 2 ) as functions of 
the intracellular concentrations C A of ATP [12,14] and C a of amino acids [13] in the 
form: 



R g i(C a ,C A ) = r g 



Ca 4 + k x 4 Ca 12 



12 



p. 12 C a s + a g 



(Al) 
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MC, C A ) = r* jf k/c J g + p n (l + (A2) 

where j s tne rate constant for autophagosome formation, with appropriate constants 
k gf p g (for ATP), y gi and ^ (for amino acids). In our simulations, the basal level of 
autophagy is suppressed or promoted by adjusting the value of r g (in units of the nor- 
mal value rg) (see Table 1). 

Autolysosome formation 

We next consider the autolysosome formation step, which consists of multiple fusions 
with lysosomes, which provide an acidic environment and a digestive function to the 
interior of the autophagosome [65,66]. In view of the experimental demonstration that 
the dynamics of autophagosome maturation depends on whether the ATP supply is on 
or off, we incorporate the intracellular ATP dependence of the step, and take the spe- 
cific rate Ri in the form [12,14]: 

where 77 denotes the rate constant for autolysosome formation, with appropriate con- 
stants I<i and pi for ATP. Note here that the possible difference between the maturation 
rates of autophagosomes from proteins/organelles Si and S 2 [51,53], which relatively 
lacks proven molecular mechanism, has not been considered in this model. Although 
extensive characterization of ATG genes has yielded critical insight into the mechanism 
of autophagy activation and its flux, so far whether and how the selective fusion of 
autophagosomes to lysosomes is controlled remain poorly elucidated [50,51,54,67-70]. 



Intralysosomal hydrolysis 

The ATP dependency of the intralysosomal hydrolysis step, which displays relatively 
low sensitivity to the ATP concentration change compared with the autophagosome 
formation step, is incorporated. In accord with experiment [12,14], the intralysosomal 
hydrolysis specific rate is taken as a function of the intracellular ATP concentration: 

*™ = n c?TkF (A4) 



with appropriate exponent 8 h and constant k^ for ATP, where is the rate constant for 
intralysosomal hydrolysis. Considering that there is little experimental evidence about 
the selective intralysosomal hydrolysis, we assume that the intralysosomal hydrolysis 
rates for autolysosomes originated from both Si and S 2 are the same. 

Protein synthesis and non-lysosomal degradation 

In addition to the autophagic process, we incorporate the (total) protein synthesis rate 
R s , depending on the amino acid concentration C a , in agreement with experiment [34]. 
Assuming that the rate grows exponentially with the intracellular ATP concentration 
C A increased to the steady-state value, we write the protein synthesis rate in the form: 
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C a exp[Q]-l (w) 

foro<c * (A5) 

with appropriate constant /c 5 for amino acid, where is the ATP concentration cor- 
responding to the maximal protein/organelle synthesis rate and r s denotes the rate con- 
stant for the protein/organelle synthesis. 

Further, non-autophagic degradation machinery such as the ubiquitin-proteasome 
system has been considered in the model. We suppose that the amount of protein deg- 
radation by autophagy constitutes up to 80% of the total amount of protein degradation 
[71]. Taking the rate of non-autophagic degradation to be 25% of autophagic degrad- 
ation, we have the rate of non-autophagic degradation (i = 1, 2): 

Rdi = \RhCu. (A6) 

where C# denotes the concentration of autolysosomes from % 
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